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We propose to use a two-species Fermi gas with the interspecies s-wave Feshbach resonance to 
realize p-wave superfluidity in two dimensions. By confining one species of fermions in a two- 
", dimensional plane immersed in the background three-dimensional Fermi sea of the other species, an 

■ attractive interaction is induced between two-dimensional fermions. We compute the pairing gap in 

' the weak-coupling regime and show that it has the symmetry of p x +ip y . Because the magnitude of 

, the pairing gap increases toward the unitarity limit, it is possible that the critical temperature for the 

Px+ipy-wave superfluidity becomes within experimental reach. The resulting system has a potential 
application to topological quantum computation using vortices with non-Abelian statistics. We 
also discuss aspects of our system in the unitarity limit as a "nonrelativistic defect conformal field 
theory (CFT)". The reduced Schrodinger algebra, operator-state correspondence, scaling dimensions 
of composite operators, and operator product expansions are investigated. 
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I. INTRODUCTION 



Experiments using ultracold atomic gases have achieved great success in realizing a new type of fcrmionic superfluids. 
i-^J , By arbitrarily varying the strength of interaction via the Feshbach resonance, the weakly-interacting BCS superfluid, 
the strongly-interacting unitary Fermi gas, and the Bose-Einstein condensate of tightly-bound molecules have been 
observed and extensively studied [1, 2], So far, the fermionic superfluids in atomic gases have been limited to s-wave 
i pairings between different spin states. Therefore the realization of p-wave superfluids in spin-polarized Fermi gases 

\ is a natural next goal in the cold atom community. In particular, a "weakly-paired" p x +ip y -wave superfluid in two 
£NJ ■ dimensions is of special interest because its vortices support zero-energy Majorana fermions and exhibit non-Abelian 
statistics [3]. As a practical application, it has been proposed to use such a system as a platform for topological 
quantum computation [4]. 

Most theoretical studies regarding the p-wave superfluids in atomic gases assume the availability of p-wave Feshbach 
resonances [4-11] (for alternative mechanisms, see Refs. [12-17]). However, experimental studies showed that the p- 
wave Feshbach molecules are unstable due to atom-molecule and molecule-molecule inelastic collisions with their 
lifetimes up to 20 ms [18-25]. This is in contrast to the long-lived s-wave Feshbach molecules where the inelastic 
collisions are suppressed due to the Pauli exclusion principle [26]. Because the decay rate of the p-wave Feshbach 
molecules is comparable to the interaction energy scale, the p-wave superfluid without additional mechanism to 
suppress the inelastic collisions will not reach its equilibrium before it decays [27, 28] . 

In this paper, we propose a novel approach to realize the p-wave superfluidity in two dimensions, without assuming 
the p-wave Feshbach resonance. The idea is to utilize a two-species Fermi gas (fermion atomic species A and B) 
with the interspecies s-wave Feshbach resonance in 2D-3D mixed dimensions [29]. Here A atoms are confined in a 
two-dimensional plane (2D) by means of a strong optical trap, while B atoms are free from the confinement and hence 
in the three-dimensional space (3D). It has been shown that the interspecies short-range interaction between A and B 
atoms is characterized by a single parameter, the effective scattering length a e ff, whose value is arbitrarily tunable by 
the interspecies s-wave Feshbach resonance [29]. The system under consideration can be set up in experiments with 
the use of the recently observed quantum degenerate Fermi- Fermi mixture of 6 Li and 40 K atoms and their interspecies 
s-wave Feshbach resonances [30, 31]. 

In such a system, we will show that the background 3D Fermi sea of B atoms induces an attractive interaction 
between A atoms in 2D. Because A atoms are identical fermions, the dominant pairing takes place in the p-wave 
channel. We will compute the pairing gap in the controllable weak-coupling regime and show that it has the symmetry 
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of Px+ipy Because the magnitude of the pairing gap increases toward the unitarity limit |a e ff| ~~ * 00 > the critical 
temperature for the p x +ip y -wave superfluidity is expected to become within experimental reach. As it is mentioned 
above, the resulting system has a potential application to topological quantum computation using vortices with 
non-Abelian statistics [3, 4]. 

This paper is organized as follows. In Sec. II, we describe the two-species Fermi gas in the 2D-3D mixed dimensions. 
In particular, we give its field-theoretical formulation in a detailed way because such a system may not be familiar to 
the cold atom community. Then in Sec. Ill, wc compute the induced interaction between two-dimensional fermions, 
the pairing gap, and its symmetry in the weak-coupling regime where we can perform the controlled perturbative 
analysis. Finally, summary and discussions are given in Sec. IV and here a very interesting analogy of the system 
investigated in this paper with the brane- world model of the universe is pointed out. Two additional materials are 
presented in Appendices. The absence of the interspecies pairing at weak coupling is shown in the Appendix A. In the 
Appendix B, we discuss aspects of our system in the unitarity limit as a nonrelativistic defect conformal field theory. 
We derive the reduced Schrodinger algebra and the operator-state correspondence in general nonrelativistic defect 
conformal field theories. We also study scaling dimensions of few-body composite operators and operator product 
expansions in our 2D-3D mixed dimensions. In particular, critical mass ratios for Efimov bound states are obtained. 



II. TWO-SPECIES FERMI GAS IN 2D-3D MIXED DIMENSIONS 



A. Field theoretical formulation 



The two-species Fermi gas in the 2D-3D mixed dimensions is described by the following action (here and below 
h = 1 and k& = 1): 



V 



2mA 



t^Aj 1pA(t,x) 

72 , Y72 



S = Jdt Jdx ip A (t, x) (id t 

+ jdt Jdx J dz^ B (t,X,z) (id t + V ^J Z 1>B(t,X,z) (1) 

+ 3o Jdt Jdx ^ A (t,x)^ B (t,x,0)ip B {t,x,0)ipA{t,x). 

Here x = (x, y) is a two-dimensional coordinate and (x, z) is a three-dimensional coordinate. ipA(t, x) is a fermionic 
field describing A atoms confined in a two-dimensional plane located at z — and tpB(t, x, z) is another fermionic field 
describing B atoms in the three-dimensional bulk space. m,A(B) is the atomic mass of A(B) atoms and the density 
of each species Ua{ B ) is controlled by the chemical potential Ha{b)- The interspecies interaction is short-ranged and 
thus occurs only on the plane at z = 0, while B atoms can propagate into the z-direction ("extra dimension") [see 
also Fig. 1]. 

go is a cutoff dependent bare coupling. Because dimensions of the fields are [iPa] = 1 and [ip B ] = § in units of 
momentum, the dimension of the coupling becomes [go] = — 1. This implies that the theory has a linear divergence as 
it is well known in the usual 3D case. However, as we will see below, the linear divergence can be renormalized into 
go and all physical observables can be expressed in terms of the physical parameter, the effective scattering length 
a e fj. We note that interactions between the same species of fermions (without the p-wave Feshbach resonance) are 
generally weak and can be neglected at low energies. 

The bare propagator of if) a field is (T ipA{t, x)ip A (t' , x'))q where the expectation value is evaluated with the non- 
interacting action. Because of the translational symmetry in the plane, its Fourier transform is given by the usual 
form: 

iG A (po,p) = - 2 -, (2) 

Po - itl + + tS 

where po is the frequency and p — (p x ,p y ) is the two-dimensional momentum. Similarly the bare propagator of ips 
field is given by (TipB{t, x, z)^ B (t', x', z'))q. We shall not perform its full Fourier transformation because once the 
interaction between if) a and ip B fields is turned on, the translational symmetry along the z-direction is lost. Instead 
it is convenient to employ the following mixed representation: 

/dp z e ip z (z-z') 
«f :, (3) 
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FIG. 1: Two-particle scattering of A and B atoms. The A atom is confined in a plane located at z — while the B atom can 
propagate into the z-direction ("extra dimension"). The interspecies short-range interaction takes plane only on the plane. 



where p z is the momentum conjugate to z — z' . We will often use the propagator where z and z' are fixed on the 
plane; z = z' = 0. In such a case, we suppress the last argument in G B (po,P] z — z') and denote it simply as 
Gb{po,p) = G B (p o ,p;0). Hereafter we shall use a shorthand notation p = (po,p). 



B. Two-particle scattering in vacuum 

We first study the two-particle scattering in vacuum (fiA = (J>b = 0) in order to relate the bare coupling go with 
the effective scattering length a e g. The scattering process of A and B atoms is schematically depicted in Fig. 1. By 
summing a geometric series of Feynman diagrams, the scattering amplitude A(p) is written as 



Mb)]-^ — - f^ tGA (p-k) l G B (k) 
igo J {2iry 

/ WE (4) 



igo J (27r) 

1 f dk 

+ i 



\ 2m. a 2m d ^ u 



We can see that the k integration is ultraviolet divergent. The usual way to regulate the integral is to introduce a 
momentum cutoff |fe| < and adjust the A^-dependence of go so that the physics does not depend on The 
integration over k leads to 

Ap) = ; 1 (5) 

i - (A* - y/z&p* - 2m ABPo - i0+) 

where M = tua + ms is the total mass and tuab = mA+m B ^ s ^ e re duced mass. By introducing the effective 
scattering length through 

1 y/m B m AB sJm B m A B ... 

o A k = — t. > y°) 

go i-K 2ira c s 

the scattering amplitude becomes cutoff-independent: 

AM = 2n 1 (7) 

^m B m AB + ^/^p 2 - 2m A BPo ~ i0+ ' 

Now the interspecies interaction is solely characterized by the effective scattering length a e g. a e g — > — corresponds 
to the weak attraction and a c s — * +0 corresponds to the strong attraction just as in the usual 3D case. |a e ff| — * oo 
corresponds to the unitarity limit where the scale-invariant interaction is achieved. In this limit, our theory (I) 
provides a novel type of nonrelativistic conformal field theories. Aspects of our system in the unitarity limit as a 
nonrelativistic conformal field theory will be elaborated in detail in the Appendix B. 

When a ff > 0, there exists a shallow two-body bound state composed of A and B atoms. Its binding energy £b, 
defined to be positive, is obtained as a pole of the scattering amplitude when the external momentum p is zero: 

^l(-£ b ,O)- 1 -0 =► Eb = - l — r . (8) 

Thus our definition of a c g in Eq. (7) coincides with that used in Ref. [29] . The two-body resonance £b — ► occurs at 
infinite effective scattering length a c g — > oo. 
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FIG. 2: Effective scattering length a e g/l as a function of the inverse bare scattering length l/a for mass ratios mA/rriB = 0.15 
(left) and mi/wiB = 6.67 (right) [29]. The vertical dotted line indicates the position of the broadest resonance. 

C. Effective versus bare scattering lengths 

The effective scattering length a c s in the 2D-3D mixed dimensions depends on the bare scattering length a in a 
free 3D space, a is arbitrarily tunable by means of the interspecies s-wave Feshbach resonance as a function of the 
magnetic field applied to the system [31]. In Ref. [29], the dependence of a e g on a was determined when the A atom 
is confined by a one-dimensional harmonic potential with the oscillator frequency u> z . Fig. 2 shows a c ff/l plotted as a 

function of l/a with I = \J being the oscillator length [29]. Here the mass ratios m^/mj = 0.15 and 6.67 are 

chosen corresponding to the physical cases of A = 6 Li, B = 40 K and A = 40 K, B = 6 Li, respectively. 

We can see that the position of the resonance in the 2D-3D mixed dimensions (|<x e ff| = oo) is shifted from the free 
space resonance (|a| = oo) to the negative bare scattering length. It is understandable that the AB bound state can 
be formed with a weaker attraction (a < 0) because of the partial confinement of the A atom. The broadest resonance 
occurs at l/a = —0.882 for niA/rnB = 0.15 and at l/a = —0.0237 for m^/rns = 6.67. In addition to the broadest 
resonance, an infinite number of confinement-induced resonances appears while they are narrower [32, 33]. 

Using one of these resonances, the effective scattering length can be tuned to any desired value — oo < a~g < oo 
by simply varying a or I. If the confinement length I is much smaller than any other length scales of the system such 
as a e ff and mean interatomic distances at finite densities, we can neglect the motion of A atoms in the confinement 
z-direction. Then the resulting system becomes the two-species Fermi gas in the 2D-3D mixture universally described 
by the action (1). 

Ref. [29] also found that the many-body system near the unitarity limit |a c ff | — > oo is stable against the formation 
of deep three-body bound states (Efimov effect) when the mass ratio is in the range 0.0351 < tiia/tub < 6.35 (see also 
the Appendix B 2). Therefore the combination of atomic species, A = 6 Li and B = 40 K (to^/tob = 0.15), can be used 
to realize the stable 2D-3D mixed Fermi gas, while the opposite combination, A = 40 K and B = 6 Li (m^/ms = 6.67), 
suffers the Efimov effect. However, because the mass ratio of the latter combination is just above the critical value, 
it may be possible that such a system becomes metastable, for example, in an optical lattice. We also note that if 
either A or B atoms are bosonic, the Efimov effect takes place for any mass ratio [29]. Thus for the stability of the 
many-body system, fermion atomic species A and B are essential. 

D. Perturbation theory at finite density 

In the limit of weak attraction a c ff — > —0, it is straightforward to develop a perturbation theory at finite densities 
(fiA> Mb > 0). The propagator of A atom is given by iGa(j>) in Eq. (2) and the propagator of B atom is given by 
iGsipiz) in Eq. (3). From Eq. (7), we find that each interaction vertex carries a small coupling constant given by 

27rm c ff 
^/m B m A B ' 

As one of applications of the perturbation theory, we compute the density distribution of B atoms in the weak- 
coupling limit a c ff — > —0. Due to the lack of translational symmetry in the z-direction, the density of B atoms is 
no longer uniform. The density of B atoms is given by ns(|z|) = (tp B (t + + , x, z)ips(t, x, z)), which is a function 
of \z\ because of the in-plane translational symmetry and the symmetry under z-parity. To the leading order in a e g, 
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FIG. 3: K(fcFs|z|) appearing in the density distribution of B atoms in Eq. (11). At z = 0, we have k(0) = l/(64.7r). 
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ns(|z|) is obtained as 

™b(M) 

, (2m B UB) 3/2 ■ 

where n B = - — g^" — i 

proportional to the effective scattering length a e ff and the density of A atoms n A ~ ( 2m ^ A ) 



dp Q dpdp z dq z 



(2tt) 3 



(9) 



(2tt) £ 



Po 



p 2 

2m b 



+ /i B + 1(5 po — 2^~r + VB + iS 



is the uniform density of B atoms in the noninteracting limit. T, B is a mean- field self-energy 

I All 



-B 



2na cS f dp dp i + 27ra G ff 



^m B m AB J (27r) 



=riA < 0. 



(10) 



^m B m AB 

We note that the dimensions of n A and are different because n A is the two-dimensional density while n B is the 

1/2 

three-dimensional density. The Fermi momentum of each species is defined through its density by kp A = (AiTn A ) 1 
and kp B = i&^ns) 1 ^ . 

The integration over po in Eq. (9) results in the following expression for the density distribution: 



n B (\z\) = n B + \a cS \k^ A k^ BX / K,(k FB \z\), 



m AB 



where n(r) is a positive function given by 



2tt 



2tt 



27r Pi - 9 Z 2 



(11) 



(12) 



k(/cfb|z|) is plotted in Fig. 3 and monotonously decreases as a function of fcps|z|. We can understand that B atoms 
are attracted to the 2D plane at z — because of their attractive interaction with A atoms confined in the plane. The 
density of B atoms away from the 2D plane approaches that in the noninteracting limit; ns(|z| — > oo) — > n B because 
the interaction is suppressed there. 



III. INDUCED INTERACTION AND p-WAVE PAIRING IN TWO DIMENSIONS 

A. Induced interaction at weak coupling 

Using the perturbation theory in the weak-coupling limit a c g — > —0, we now determine the interaction between two 
A atoms in 2D induced by the existence of the 3D Fermi sea of B atoms. Because we are interested in the intra-species 
pairing of A atoms, we consider their back-to-back scattering. To the leading order in a e s, the induced interaction 
between A atoms V; n d is described by the Feynman diagram depicted in Fig. 4 [16], which is written as 



1 . . If -2wia e g 

2 2 \y/m B m AB 



dkadk 



iG B {k+p- q)iG B {k). 



(13) 
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FIG. 5: v(w) appearing in the induced interaction in Eq. (15). 



The integration over ko leads to 



Mnd(p, q) 



(27ra e ff) f dkdk z dk' 



m B m A B 



q( (fc+p-q) 
\ 2ms 



(2^) 4 



- Ms) # Mb ~ 



2m B 



-) ^(mb 



(fc+p-q) 2 + fcj 
2m B 



fc 2 +fc; 2 
2m B 



MB 



(fc+p-q) 2 +fcj 

2m_B 



k 2 + k'? 
2m B 



(k+p-q) 2 +k? 
2m B 



+ Po - 9o - «0 H 



(14) 



For the gap equation at weak coupling, we will need the induced interaction in which both incoming and outgoing 
momenta are on the 2D Fermi surface \p\ = \q\ = fcpAj and hence, po = 9o = 0- I n such a static limit, we can perform 
the remaining integrations analytically and obtain 



V ind (p, q) 



27ra c 2 ff fcp B 
mAB 



2A- 



FB 



where u(u') is a continuous function given by 



(15) 



v(w) = 



Vw 2 - 1 + (2- w 2 ) 



arcsm w 



4tt 



to < 1 



w > 1. 



(16) 



The nonanalyticity of v(iu) at it) = 1 is due to the sharp Fermi surface of -B atoms. The function v(w) is plotted in 
Fig. 5 and is negative everywhere indicating that the induced interaction between A atoms is attractive. Thus an 
intra-species pairing in the two-dimensional plane is expected to occur. Because A atoms are identical fermions, the 
dominant pairing takes place in the p-wave channel, as we will see below. 

At this point, we should point out that the interspecies pairing between A and B atoms is unlikely in our system 
(except deep in the BEC regime a e g — > +0) because they live in different spatial dimensions. B atoms can always 
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escape from the 2D plane in which A atoms are confined into the z-direction ("extra dimension") and there the 
interspecies interaction is turned off. Actually, as we will show in the Appendix A, the absence of the interspecies 
pairing can be confirmed at weak coupling a c g — > — 0. Hereafter B atoms are treated as a background to induce the 
attraction between A atoms and we investigate the intra-species pairing of A atoms in 2D. 



B. Gap equation 

Once the induced interaction between A atoms Vind(p, q) is obtained, the pairing of A atoms in 2D is described by 
the BCS-type Hamiltonian: 



HA = /(^fe-^)^' [(P 



(17) 

+ 1 / ^i(H a(* - «) y^MH **(f +*) • 

where 4>a{p) is the Fourier transform of %Pa(x). We note the property Vind(p, q) = Vj n d(g,J>). The pairing gap of A 
atoms A p is defined to be 



(2TrfS(k)A p = JJS_ Vind (p,q) 



M-z-qiM ( 18 ) 



Because of the Fermi statistics of A atoms, the pairing gap has to have an odd parity; A_ p = — A p . The standard 
mean-field calculation leads to the following self-consistent gap equation: 

Ap = / (Sp Mnd(P ' Q) 2l| [1 " 2nF{E " )] ■ (19) 



Here E p = y (2^7 — ' u - 4 ) I^pI 2 * s tnc q uas ip ar ticlc energy and n-p(E p ) = 1/ (e Ep / T + l) is the Fcrmi-Dirac 

distribution function at temperature T. 

The gap equation (19) is a nonlinear integral equation in terms of the pairing gap A p . However, it becomes a 
linear integral equation near the critical temperature T — > T c because one can set A p — > in E p . In such a case, 
A p = e M *AW with an odd integer I being the orbital angular momentum solves the gap equation and the critical 
temperature T c is determined by the equation 



-tanh — m . (20) 



Here A £ is an energy cutoff and TV^ = ^ is the density of states of A atoms at the Fermi surface. — V^ d ^ is 
the partial-wave projection of the induced interaction given by 



(0 = 2na^ B rtB_ J ll - cosfl \ 2na^ B (1) / W 

md m AB Jo 7T 1 ^ ^fcpsV 2 y ITlAB \k FB/ ' 

When the projected interaction is attractive A/^V^j < 0, it is easy to solve Eq. (20) and we find 



(naV") 



T d) 

-^-~exp( — — ). (22) 

k 2 

where the energy cutoff is chosen to be the order of the Fermi energy of A atoms; A e ~ sfa = 2mX' Because one can 
confirm that the induced attraction (21) is strongest in the p-wave channel, we have the highest critical temperature 
for |/| = 1; T c (1) > T c (|Z| - 3) . Therefore we can neglect the coupling between different partial waves and concentrate on 
p-wave pairings. 
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C. Pairing gap and its symmetry 

We now solve the gap equation for the p-wave pairing at zero temperature T = 0. We parameterize the angle 
dependence of the pairing gap as A p = f(6p)Af, where p — fcFyt(cos# p , sin# p ) and f(9p) = b + e l0 p + b_er l0 v with 
\b + \ 2 + |6_| 2 = 1. For example, b + = 1 and = corresponds to a p x +ip y -wsve pairing and b + = fo_ = l/\/2 
corresponds to a p^-wave pairing. Substituting A p = f(9p)Af into the gap equation (19) at T = 0, we obtain 



o Jo 



2tt V^ + I/^A/P (23) 



/ 27r'^*" VI/(^)A/ 
Thus we find that the modulus of the pairing gap is given by 

T^expf-i-^- rg|/(0)| 2 ln|/(0)|) . (24) 

The angle dependence of the pairing gap f(9) is determined so that the ground state energy is minimized [34]. 
Because the gain of energy density due to the condensation is given by 

(^-(^l^^-^IA,! 2 , (25) 

the ground state energy is minimized when |A/| is maximized. From Eq. (24), we can show that the maximum |A/| 
is achieved when f{9) = e ±l0 corresponding to the p x ±ip y -wave pairing. Therefore the pairing gap realized in our 
system becomes 

^-^"'expf^—). (26) 

We note that the pairing symmetry p x ±ip y is favored because of the isotropy of the induced interaction; V^} = V^ 1 ^ . 
This is in contrast to the p-wave Feshbach resonance where the interatomic interaction can be anisotropic due to the 
magnetic dipole-dipole interaction [35]. In such a case, some parameter-tunings are necessary to realize the p x ±ip y - 
wave pairing [8, 10]. 



D. Optimizing the pairing gap 



In order for the experimental realization of the proposed p x +ip y -wave superfluidity, the critical temperature (22) 
and the magnitude of the pairing gap (26) have to be large enough. Thus we look for a condition in which ~ 
I A p | ~ £ya exp 1/ {NaV-^J^ is maximized. From Eq. (21), NaV^J is given by 



mAB \k F Bj 



(27) 



with the negative function i/ '(fepAAFB) plotted in Fig. 6. Because the induced interaction in the p-wave channel is 
attractive NaV^J < 0, one would like to minimize NaV^J. 

One possible way to optimize the pairing gap is to control the densities of A and B atoms [16]. Because our 
perturbative calculation relies on the smallness of |a e fffcFs|, we fix a c ffkpB and vary the ratio in the two Fermi 
momenta k~FA / ^fb • We find that the function (feEA/ftpe) has an minimum = —0.0452 at /cea / &fb = 1-75 (see 
Fig. 6). Thus, within our perturbative calculation, the maximum pairing gap becomes 

A™ ax ,., / 22.1 m AB \ 

- e±^exp -— -** . (28) 



£fa \ m A a^ s kj B 



9 




The pairing gap can be further enhanced by changing the mass ratio m^/m^. Because of the factor vhab/^a in 
the exponent, the larger mass of A atoms in 2D increases the pairing gap. For example, the combination of atomic 
species A = 40 K and B — 6 Li has ra^/mj = 6.67, and hence, the pairing gap becomes 

A Max / 3 2Q \ 

e ±% ex P (-45r). (29) 



cffTB 



Now the exponential factor is not hopelessly small. If one could extrapolate our perturbative result to |a e fffcEs| ~ lj 
we would have expl IttV ) ~ 0.04. Furthermore, in the unitarity limit |a ff I — ► oo, the pairing gap is expected to 

V "eff FB J 

be the same order as the Fermi energy; A p ~ e^^e-pA- Therefore it is possible that the critical temperature for the 
p x +ipy-~w&ve superfluidity becomes within experimental reach, in particular, near the unitarity limit. 



E. Nonperturbative approaches near the unitarity limit 

So far, we have performed the controlled perturbative analysis in the weak-coupling regime a c g — ► —0. An important 
quantitative question is how high the critical temperature 7c can be near the unitarity limit |a e fj| — > oo. Ideally 
one would like to answer this question by employing quantum Monte Carlo simulations while they will suffer fermion 
sign problems because of the intrinsic asymmetry between A and B atoms in our 2D-3D mixture. Instead it is 
possible to estimate Tc by using nonperturbative analytical methods such as the e expansion [36-39] and the l/N 
expansion [40, 41]. 

The application of the l/N expansion technique to our 2D-3D mixture is straightforward. We generalize the two- 
species Fermi gas in Eq. (1) to a (2A^)-species Fermi gas in which N species live in 2D while the other N species live 
in 3D. The interaction among them occurs on the 2D plane and is assumed to be the Sp(27V)-symmetric form [40, 41]. 
Then we utilize the small parameter l/N <C 1 to perform systematic expansions. 

Here we comment on the application of the e expansion technique to mixed-dimensional systems. Suppose A and B 
atoms live in cIa- and ^-dimensional spaces, respectively, where the former space is a subset of the latter space with 
Aa < d>B- Such a system is described by the action analogous to Eq. (1). Now the dimensions of the fields change to 
[i/)a] — <1a/2 and [V>b] = rfs/2 in units of momentum, and thus, the dimension of the coupling becomes [go] = 2 — ds- 
As far as [go] > —2 is satisfied, the theory is renormalizable [37]. We note that [go] depends only on the bulk spatial 
dimension ds indicating that ds plays a central role in the e expansion. In the general combination of the spatial 
dimensions, one can study the two-particle scattering in vacuum as it was done in Sec. II B. Using the dimensional 
regularization, the scattering amplitude A(p) at the scale-invariant unitarity point is found to be 

(jHB_Y A/2 (2n_\ dB/2 

«P) 1 7, [mB) ( 3 °) 

where p is the o^-dimensional momentum. We can see that the scattering amplitude vanishes in the limits of ds ^4 
and ds — * 2 indicating that these two spacial dimensions correspond to noninteracting limits. Accordingly we can 
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FIG. 7: Proposed phase diagram of a two-species Fermi gas in the 2D-3D mixture as a function of the inverse effective scattering 
length (a £ffepA) _1 at zero temperature. There is a quantum phase transition at (a e ff A;fa) _1 — 0(1) from the 2D-3D mixed 
Fermi gas with the 2D p^+ip^-wave superfluidity (2D p-SF) to the 2D Bose-Einstein condensation of localized s-wave molecules 
(2D BEC). 



develop systematic expansions in terms of e = 4 — <ie <§; 1 and e = Cl around those special dimensions [36-39]. 

We note that the other spatial dimension d^(< ds) is arbitrary in this approach. 

The estimation of near the unitarity limit using the above nonperturbative approaches will be left for future 
works. 



IV. SUMMARY AND DISCUSSIONS 



In this paper, we presented theoretical prospects to realize the p-wave superfluidity in two dimensions by using 
a two-species Fermi gas (fermion atomic species A and B) with the interspecies s-wave Feshbach resonance. By 
confining A atoms in a 2D plane immersed in the background 3D Fermi sea of B atoms, an attractive interaction is 
induced between A atoms. Because A atoms are identical fermions, the dominant pairing takes place in the p-wave 
channel. In the weak-coupling regime a c g — > —0 where the controlled perturbative analysis is available in terms of the 
effective scattering length, we computed the pairing gap and showed that it has the symmetry of p x +ip y . Because 
the magnitude of the pairing gap increases toward the unitarity limit |a e ff| - > 00 , the critical temperature for the 
Px+zpy-wave superfluidity is expected to become within experimental reach. As it is mentioned in the Introduction, 
the resulting system has a potential application to topological quantum computation using vortices with non-Abelian 
statistics [3, 4]. 

It is worthwhile to clarify what happens deep in the BEC regime a c g — > +0 in our 2D-3D mixture. In this limit, 
A atoms in 2D capture B atoms to form tightly-bound molecules and the resulting system consists of the molecules 
localized on the 2D plane plus excess A or B atoms. When the size of the molecules ~ a c g becomes smaller than the 
mean interatomic distance in 2D ~ k^ A , the molecules behave as two-dimensional bosons and therefore the ground 
state will be a 2D Bose-Einstein condensate of the s-wave molecules. Consequently, there has to be a quantum phase 
transition from the 2D-3D mixed Fermi gas with the 2D p x +ip y -w&ve pairing [— oo < (aefffovi) -1 ^ 0(1)] to the 
2D Bose-Einstein condensation of the s-wave molecules [0(1) < (a e fffcEA) _1 < +oo]. The proposed phase diagram 
as a function of the inverse effective scattering length is shown in Fig. 7. These two phases can be distinguished by 
radio-frequency spectroscopy experiments. In the 2D-3D mixed Fermi gas with the 2D p x +ip y -w&ve pairing, A atoms 
are fully gapped while B atoms remain gapless. On the other hand, in the 2D Bose-Einstein condensation of the 
s-wave molecules, both A and B atoms are fully gapped. 

Readers may wonder why we did not consider a system in which both A and B atoms are confined in a two- 
dimensional plane to realize the induced p-wave superfluidity in 2D. In this case, A and B atoms always form bound 
molecules in 2D and thus the ground state of the system tends to be an s-wave paired state. In order to break the 
interspecies s-wave pairing, one needs to weaken the interspecies attraction with a large density imbalance intro- 
duced [42] . This would be a disadvantage in order to achieve a high critical temperature for the p-wave superfluidity 
in 2D. 

A remarkable aspect of our two-species fermions in the 2D-3D mixed dimensions is that the system in the unitarity 
limit |a c ff| ^ oo is described by a nonrelativistic defect conformal field theory, which is a novel class of quantum field 
theories that has not been paid attention to so far. We elaborated this aspect in detail in the Appendix B. 

Finally, it is very interesting to point out the analogy of the system investigated in this paper with the brane- 
world model of the universe. In the brane-world scenario, the ordinary matter is considered to be confined in 
a three-dimensional space (brane) embedded in higher dimensions (bulk) where gravitons can propagate [43]. The 
gravitational force between matters is induced by the exchange of the graviton. Similarly, in our system, the interaction 
between A atoms confined in the 2D plane ( "2D brane" ) is induced by the exchange of B atoms in higher dimensions 
("3D bulk"). Within this fascinating analogy, our two-species Fermi gas in the 2D-3D mixed dimensions can be 
regarded as a brane world in cold atoms! 
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APPENDIX A: ABSENCE OF INTERSPECIES PAIRING AT WEAK COUPLING 



Here we show the absence of the interspecies pairing between A and B atoms in the 2D-3D mixed dimensions in 
the weak-coupling regime a c g — > —0. We consider the Nambu-Gor'kov-type propagator in the 2x2 matrix form: 



ig(t-t',x-x') 



( (TipA(t, x)ip A (t' , x')) (T^ A (t,x)^ B (t',x',0)) 
\(T 4(t, x, 0)^(f, *')> (T^ B (t, x, 0)^ B (f, x', 0)) 



(Al) 



In the mean-field approximation, the above propagator in the momentum space becomes 

Gip) 



, G-\-p) 

_ I G-/(p)G B \-p) + \M 2 G- A \p)G B '(-p) + \<t>o\ 2 



(A2) 



,G- A 1 (p)G B 1 {-p) + \M 2 G- A 1 {p)G B 1 (-p) + \M 2 l 

where 0o = (go V'sC*, x, 0)Vu(i, x)) is a condensate determined by the self-consistent gap equation: 



<h _ _. f dm 

go ~ J (2tt 



dp ~ 

? 



m B 



dpodp 



(A3) 



(2tt)3 



(ipo 



+ /m) \Jipo + — 



2m a 



In the last line, we analytically continued po to the imaginary frequency; po 
length via Eq. (6), we obtain the following renormalized gap equation: 



2rn B 

ipo. Introducing the effective scattering 



y/2m A B _ f dppdp 
~ _ J ~J^f 



2ira, 



2m B 



Ms 



(A4) 



For simplicity, we shall consider the equal masses uia = tub = m and equal chemical potentials fiA = Ms = M 
where the interspecies pairing is guaranteed in the usual 3D case. However, in the 2D-3D mixture, we can see that 

2 

the right-hand side of Eq. (A4) does not have any singularity around the Fermi surface, po ~ and <~ fi, in the 
limit |0o | —> 0, and hence, the integral is bounded from above. This can be understood as an absence of the Cooper 
instability because of the intrinsic "mismatch" between the 2D and 3D Fermi surfaces. Therefore in the weak-coupling 
regime a e s — > — 0, the gap equation (A4) does not have a nontrivial solution showing that there is no interspecies 
pairing between A and B atoms. 



APPENDIX B: ASPECTS AS A NONRELATIVISTIC DEFECT CONFORMAL FIELD THEORY 

As we mentioned in Sec. II B, two-species fermions in the 2D-3D mixed dimensions in the unitarity limit |a e ff| — * oo 
(at zero density and zero temperature) provide a novel type of nonrelativistic conformal field theories (CFTs) 1 . In 



1 As far as we know, there arc three basic ingredients to construct interacting nonrelativistic CFTs; l/J? 2 -type interactions, zero-range 
interactions at resonance [44], and interactions due to fractional statistics in two dimensions [45]. Combinations of these interactions 
also work [46, 47]. In addition to those field-theoretical constructions, gravity dual descriptions of different classes of nonrelativistic 
CFTs have been recently proposed [48, 49]. It would be interesting to investigate gravity duals for nonrelativistic defect CFTs imitating 
the situation studied in this paper. 
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TABLE I: Summary of eight classes of nonrelativistic (defect) CFTs with zero-range and few-body resonant interactions 
proposed in Ref. [29]. Two-species fermions in pure 3D are the well-known case of nonrelativistic CFT [44, 52]. In all cases 
below, the coupling of zero-range and few-body interaction term has the dimension [go] = — 1 and is tuned to the resonance. 



Nonrelativistic (defect) CFT 



Spatial configurations 



Symmetries other than H, D, C, M. 



2 species in pure 3D 

2 species in 2D-3D mixture 

2 species in 1D-3D mixture 

2 species in 2D-2D mixture 

2 species in 1D-2D mixture 

3 species in 1D-1D-1D mixture 

3 species in 1D 2 -2D mixture 

4 species in pure ID 



xa = xb — (x,y,x) 

x A = (x,y) x B = (x,y,z) 

x A = (2) x B = (x,y,z) 

x A = (x,z) x B = {y,z) 

x A = (2) x B = (x, y) 

x A = (x) x B = (y) x c = (2) 

x A = x B = (x) xc = {x,y) 

xa = xb = xc = xd = (x) 



Pi, Ki, Jij with i,j = x,y,z 
Pi, Ki, Jij with i,j = x,y 
P K 7 

1 z , J^z , *Jxy 

P„ K z 

Jxy 

None 
P K 

1 x, l^x 

P K 

1 x, l^x 



this system, the three-dimensional translational, rotational, and Galilean symmetries in the bulk space are broken to 
the two-dimensional symmetries while scale and conformal invariance are preserved. Regarding the two-dimensional 
plane as a deject in the three-dimensional bulk space, our system can be thought of a nonrelativistic counterpart of 
defect/boundary CFTs [50, 51]. In Ref. [29], more classes of nonrelativistic defect CFTs with zero-range and few-body 
resonant interactions have been proposed and are summarized in Table I. Here we discuss aspects of our system as 
a nonrelativistic defect CFT (abbreviated as NRdCFT). First we derive the reduced Schrodinger algebra and the 
operator-state correspondence in general nonrelativistic defect CFTs. Then we study scaling dimensions of few-body 
composite operators and operator product expansions in our 2D-3D mixture. In particular, the critical mass ratios 
for the Efimov effect are obtained. 



1. Reduced Schrodinger algebra and operator-state correspondence 



Here we derive the reduced Schrodinger algebra and the operator-state correspondence in general nonrelativistic 
defect CFTs. For defmiteness, we consider systems with two species of particles because the generalization to more 
species is straightforward. Define the mass densities 



and the momentum densities 



m A (x A ) = m A ip A (x A )tp A (x A ) 
m B (x B ) = m B iP b {xb)iPb{xb) 

Za{x a ) = -^ip A (x A )V A ijj A (x A ) 
Zb{xb) = -^ B {x B )V B Tp B (x B ). 



(Bl) 



(B2) 



Here x A (Va) is a d^-dimcnsional coordinate (derivative) and x B (Vb) is a ds-dimensional coordinate (derivative). 
We assume that the intersection of the spaces in which A and B particles live exists and includes the origin x A = x B = 
0. For example, in our 2D-3D mixture, we have x A = (x,y) and x B — (x,y,z), while in general the (iyt-dimensional 
space may not be the subset of the (ie-dimensional space such as in the 2D-2D and 1D-2D mixtures in Table I. We 
suppress the argument of time when we denote the operators ip A (t, x) and ip B (t, x, z) at t = 0. 

We consider commutation relations of the following set of operators in general mixed dimensions: the Hamiltonian 



H = dx A 



V A tp\(x A ) ■ V A ip A {x A ) 



+ / dxi 



V B ip B {x B ) ■ V B tp B (x B ) 



2m A J 2m B 

dx A I dx B xp A (x A )ipl,(x B )V(x A ,x B )xp B (x B )i; A (x A ), 



(B3) 



the dilatation operator 



D = j dx A x A ■ % A (x A ) + j dx B x B ■ Z B {x B ), 



(B4) 
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TABLE II: Full Schrodinger algebra in d spatial dimensions Sch(d) taken from Ref. [52]. The values of [X, Y] are shown 
below. The commutators of M and Jij with other operators are given by [M, any] = [Jij, D] = [Jij, C] = [Jij, H] = 0, 
[Jij, Jki] — i{5ikJji + SjiJik — SaJjk — SjkJu), [Jij, Pk] = i(SikPj — SjkPi), and [Jij, Kk] = i{&ikKj — SjkKi) with i,j = l,...,d. 



X\Y 




Kj 


D 


C 


H 


Pi 





—iSijM 


~iP l 


-iKi 





Ki 


iSijM 





iKi 





iPr 


D 


iPj 


-iKj 





-2iC 


2iH 


C 


iKj 





2iC 





iD 


H 





-iPj 


-2iH 


-iD 






and the special conformal operator 

c =\ J dx A x\m A {x A ) + i Jdx B Xgm B (x B ). (B5) 

D and C are the generators of scale transformation x A (b) — * eXx A(B), t — * e 2A i an d conformal transformation 
^AfB) - * + Af), i — » t/(l + At), respectively. The commutation relation 

[D, C] = -2iC (B6) 

can be checked by a direct calculation. By using the continuity equation [H, mA(xA)} = i^A ■ Z A {x A ) and the same 
with A — > B, we can show 

[if, C] = -i£>. (B7) 

Finally, if the interparticle interaction V(x A ,x B ) is scale invariant (for example, \x A — a3s| _2 -type interactions or 
zero-range and infinite effective scattering length interactions proposed in Ref. [29]), we obtain 

[D, H] = 2iH. (B8) 

If the system has unbroken translational, rotational, and Galilean symmetries, the corresponding generators, namely, 
the momentum operators 

Pi = j ' dx A Z Ai (x A ) + j dx B ZBi{x B ), (B9) 

the angular momentum operators 

= Jdx A [ XAiM XA) ~ x aM xa)] + Jdx B [ XB a Bj ( XB ) - x bM x B )] , (BIO) 
and the Galilean boost operators 

Ki = J dx A x Ai m A {x A ) + Jdx B x Bi m B (x B ), (Bll) 

together with the above H, D, C, and the mass operator 

M = J dx A m A (x A ) + J dx B m B (x B ) (B12) 

form the (reduced) Schrodinger algebra [52] (see Table II). Various classes of the reduced Schrodinger algebra are 
possible depending on spatial configurations of defects as shown in Table I. For example, in our 2D-3D mixture, there 
are planer translational, rotational, and Galilean symmetries preserving the location of the 2D defect at z = and 
hence we can take Pi, Ki, and Jij with i,j — x, y. In some cases such as the 1D-1D-1D mixture with three species of 
particles, all translational, rotational, and Galilean symmetries are broken by defects and thus only H , D, C, and M. 
form the reduced Schrodinger algebra. We note that the symmetry transformations in the 2D-3D mixed dimensions 
are not equivalent to those in the usual two dimensions although they have the same Schrodinger algebra Sch(d = 2). 
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This is because the scale and conformal transformations generated by D and C in Eqs. (B4) and (B5) involve the 
^-direction perpendicular to the 2D defect. 

It is useful to introduce a notion of primary operators. Consider a local operator O(xab) composed of iPa(xa) and 
*Pb(xb) operators where xab = xa = xb is a coordinate on the intersection of the <1a- and e^-dimensional spaces 
including the origin. 0{xab) is also called a defect operator because it lives on the defect. The local operator O is 
said to have a scaling dimension Ao and a mass Mo if it satisfies 

[D, O(0)} = iA o O(0) and [M, O(0)] = M o O(0). (B13) 

Furthermore when O at origin commutes with C and Kj (if Kj exists), 

[C, 0(0)} = [Kj, 0(0)} = 0, (B14) 

such a operator is called a primary operator. Starting with the primary operator 0(xab), one can build up a 
tower of local operators by repeatedly taking its commutators with H and Pj (if Pj exists) [52]. For example, 
[H, 0(xab)} = —id t O(xAB) is a local operator having the scaling dimension Ao +2 and [Pj, 0(xab)} = idjO(xAB) 
is a local operator having the scaling dimension Aq + 1 . 

We are now ready to show the operator-state correspondence in nonrelativistic defect CFTs. Consider the state 

|*o) =e- H / w O t (0)|0), (B15) 

where O is a primary operator. Then it is easy to show that [^o) is an energy eigenstate of the oscillator Hamiltonian 
H osc = H + oj 2 C with an energy eigenvalue Aqlo: 

#osc|*o) = {H + uj 2 C) e- H/u Ol\0) = e- H/bJ {u?C -iuD)O^\0) =uA \^o)- (B16) 

We note that the external potential term ui 2 C in H osc represents a e^^-dimensional harmonic potential for A(B) 
particles with (La and rig being different spatial dimensions in general. By further acting a raising operator 

L^ = — -uC + iD (B17) 

to the primary state l^o), we can generate a semi-infinite ladder of energy eigenstates (L^^^o) with n = 
0,1,2,... [52]. Their energy eigenvalues are given by (Ao + 2n)uj and can be interpreted as excitations in the 
breathing mode [53]. If Pj and Kj exist, one can make another raising operator 

«5 = i + Vl^ < B18 > 



which generates energy eigenstates (Q^) n \^o) with energy eigenvalues given by (Ao +n)oj. They correspond to 
excitations in the ccntcr-of-mass motion. The lowering operators L = ^ — loC — iD and Qj = -^== — i^f^Kj 
annihilate the primary state; L\^>a) — and Qj\^>o) = 0. 

Generalizations of other properties discussed in Ref. [52] also hold in our nonrelativistic defect CFTs. In particular, 
the two-point correlation function of the primary operator O is determined up to an overall constant in terms of its 
scaling dimension Aq and its mass Mo [52, 54]: 



(TO(t,x AB )O\0,0)) ozt- A °cxp^-iM \-^4^j . (B19) 



2. Composite operators and anomalous dimensions 

We now turn to our specific nonrelativistic defect CFT, namely, two-species fermions in the 2D-3D mixed dimensions 
(1) in the unitarity limit |a c ff | — » oo. Here we study various primary operators and determine their scaling dimensions. 
The simplest primary operators are one-body operators ^Pa(x) and ^Pb(x,0) whose scaling dimensions are trivially 
A^ A — 1 and A^ B =3/2, respectively. 

A nontrivial primary operator is the two-body composite operator 

0(x) = lim \y - x\i) B (y, 0)^ A {x). (B20) 
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FIG. 8: Feynman diagrams to renormalize the three-body composite operators <j)ipA(B)- The solid lines are the propagators of 
tpA and ips fields and the dotted lines are the propagators of <j> field. The shaded bulbs represent the vertex function Z A (b){p)- 



The presence of the prefactor \y — x\ guarantees that matrix elements of the operator <j>(x) between two states in the 
Hilbert space are finite. Thus its scaling dimension becomes 

A* = A^ A +A* fl -l = |. (B21) 

This result can be conformed by computing the two-point correlation function of (f> and comparing it with Eq. (B19): 

(T 0(4,*)^ ( Oj 0)) = -i j^.^-ivot A{p) a t -s/2 ex J m \^\ (B22) 

Here A{p) is the two-particle scattering amplitude given in Eq. (7) with |a e ff| — > oo. The <j) field can be also interpreted 
as an auxiliary field that appears when we decompose the four-Fermi interaction term in the action (1) using the 
Hubbard- Stratonovich transformation; <f>(x) = godi B (x,0)ipA(x). For the later use, we denote the Fourier transform 
of the above <f) propagator as iD(p) = —iA{p). 



a. AAB three-body operators 



We then consider three-body composite operators. A three-body operator composed of two A atoms and one B 
atom with zero orbital angular momentum I = is 

O l £°l(x) = Z^H^Aix), (B23) 

where is a cutoff-dependent renormalization factor. We study the renormalization of the composite operator (fnpA 
by evaluating its matrix element (Q\<jnl)A{x)\p, —p). Feynman diagrams to renormalize 4>i/ja is depicted in Fig. 8. The 
vertex function Z A (po,p) m Fig. 8 satisfies the following integral equation: 



Z A ( Pa ,p) = 1 - i J G A {-k)G B {k + p)D{k)Z A {k , k) 



dk 1 1 f k 2 \ (B24) 

- Z A\ —z i fe 



mAB J ^ /(|±£li + fc^_p _ ? o+ + V 2m A 

y 2m B 2m a ^ u V 2M 2m a 

where we used the analyticity of Z A (ko, k) on the upper half plane of feo- The minus sign in front of the second term 
comes from the Fermi statistics of A atoms. When we set po = — 2mX' ZA ^P) = %a(— >pj satisfies 

ZA(p) = 1 _^A [ f 1 1 



Because of the scale invariance and in-plane rotational symmetry of the system, we can assume the form of z A (p) 
to be 

za(p)= Kt) 7+1 ' (b26) 
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where A is a momentum cutoff and \ is an unknown constant. The renormalization factor becomes oc A 1 and 
thus 7 + 1 is the anomalous dimension of the composite operator (jripA ■ (Here 7 and 7; appearing later are defined so 
that they coincide with the definition of the scaling exponents used in Ref. [29].) Once 7 is determined, the scaling 

dimension of the renormalized composite operator O aa ° b is given by 



= A + + 7 + 1 = I + 7. ( [327 ) 



Substituting the form (B26) into Eq. (B25), we obtain 



|PP +1 - XA- 1 ^K=^ Ag^ ; , ffi = = ■ (B28) 



m AB + 1 •> ( 2?r ) 2 v /^-(fc + p) 2 + fe 2 +p 2 



In order for the integral to be infrared finite, Re (7) > —2 is necessary. Also, in order to be able to take the limit 
A — > 00, Re (7) < 1 is required. In the infinite cutoff limit A — > 00, 7 satisfies the following equation: 

1 = ^ r° dfc r ^ j (B29) 

\m B J m B \m B J 

where k = |fe|/|p| and cos 6* = k ■ p. Here the integral is understood to be evaluated where it is convergent — 2 < 
Re (7) < — 1 and then analytically continued to an arbitrary value of 7. 

Similarly, for general orbital angular momentum I, we consider the following three-body composite operator: 

1 

0^ AB (x) = Zl 1 Y J ^{d x + ^ y ) J ^(x)(d x +id y ) l ^^ A (x). (B30) 
3=0 

In order for 0^ AAB to be a primary operator ([-fQ, = 0), the coefficients Cj have to be chosen so that 

E * ( P + — g) j (~P+ ^ q) ^ « P l (B31) 
j^y V M + m A J V M + m A J F 1 ; 

being independent of the momentum q conjugate to the center-of-mass motion. In the important case of I — 1, we 
easily find ci = — ^j^cq. If we denote the anomalous dimension of such a composite operator as 7; + 1 — I, it is 
straightforward to show that 7/ satisfies 

m A r d Q cos(^+i (B32) 



m ^VW+lJo J° * J(n* + i) k 2 + ^kcos9+(^ + l) 

y \m B J m B \m B ) 



The integration over k leads to the result shown in Ref. [29]. The scaling dimension of the renormalized composite 
operator O aab is given by 

A^ AB = A^ + A^ + I + ( 7; + 1 - /J = I + lh (B33) 

The anomalous dimensions 7; obtained by solving Eq. (B32) in the s-wave channel I = and the p-wave channel 
I = 1 are plotted in Fig. 9 as functions of the mass ratio m A /mB- For I = 0, 70 increases as m A /mB is increased 
indicating the stronger effective repulsion in the s-wave channel. On the other hand, for I = 1, 71 decreases with 
increasing m A /rriB and eventually becomes complex as 71 = — | ± ilm(7i) when m A /mB > 6.35111 [29]. (For 
comparison, the Born-Oppenhcimcr approximation predicts the critical mass ratio to be m A /m,B ~ 6.21791.) In this 
case, using the scaling dimension A^^ = 2 ±21111(71) and Eq. (B19), the two-point correlation function of O aab is 
found to behave as 



J dtdxe- tp - x+vpat {TO{t,x)0\Q,Q)) = E 6 ± 



P 2 

-Pa- i0+ 



Am A + 2m,B 



oc sin 



±i Im(7i) 

(B34) 



, , / P 2 ~ (4m A + 2m B )po - i0 + , 
Im( 7 i) In ( — ) + Lp 
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FIG. 9: Anomalous dimensions ji for O aab in the s-wave channel I = (left) and the p-wave channel 1 = 1 (right) as functions 
of the mass ratio m A /mB [29]. In the right panel, the real part of 71 (solid curve) and its imaginary part with a negative sign 
(dashed curve) are plotted. 



Now the full scale invariance is broken to a discrete scaling symmetry, 

p -> e VIm(7i) p and pQ _> e 27r/Im(7i) pQj ( B35 ) 

which is a characteristic of a renormalization-group limit cycle [55]. This implies the existence of an infinite set of 
discrete bound states in the p-wave AAB three-body system. The energy eigenvalues form a geometric spectrum as 
E n+1 /E n = e -2^/|Im(7i)l an d tne y are known as Efimov bound states in the usual 3D case [56]. Because the system 
develops deep three-body bound states, the corresponding many-body system cannot be stable toward collapse. 

We note that an interesting thing becomes possible in the range of mass ratio 2.32780 < m^/ms < 6.35111 [29]. 
Here the anomalous dimension is — § < 71 < —\ (see the right panel in Fig. 9) and thus the scaling dimension of the 

three-body composite operator 0^7^ becomes 2 < < 3. Therefore a new three-body interaction term 

%body =9i J dtdx 0t (t, x)0{t, x) (B36) 

becomes renormalizable because now the coupling has the dimension —2 < \g{\ < [47]. The action (1) with S3-body 
added defines a new renormalizable theory. In particular, when g\ is tuned to an AAB three-body resonance, the 
resulting system provides a novel nonrelativistic defect CFT describing two-species fermions with both two-body (AB) 
and three-body (AAB) resonances in the 2D-3D mixture [29]. 



b. ABB three-body operators 

A three-body operator composed of one A atom and two B atoms with zero orbital angular momentum I = is 

0%^{x) = Zl 1 ^{x)^ B {x 1 Q), (B37) 

where Z\ is a cutoff-dependent renormalization factor. We can study the renormalization of the composite operator 
(jyipB by evaluating its matrix element (0\(fnpB(x)\p, —p)- Feynman diagrams to renormalize 4>i/jb is depicted in Fig. 8. 
The vertex function Zb(po,p) in Fig- 8 satisfies the following integral equation: 



Zb(poiP) = l-i J G B {-k)G A {k + p)D{k)Z B {k ,k) 



. 2 fdkdk z 1 1 „ ( k 2 + k? \ ( B38 ) 
1 \ / / 7Z so , s, rr—ro , Z B , k 



m AB V m B J (2tt) 3 ih+E^L + _ po _ i0 + Jy fc^+fcj \ 2m i 



2m A 2m B ™ \J jM t 

where we used the analyticity of Zsiko, k) on the upper half plane of fco- The minus sign in front of the second term 

2 . 2 / 2 1 , 2 \ 

comes from the Fermi statistics of B atoms. When we set po = ~ P 2m ^ 1 z b{p,Pz) = Zb( — P 2m ^ iP) satisfies 

AirrriB f dkdk z 1 1 

Z B (p, Pz ) = 1 ~— ] (2?r) 3 ^ {k+p) 2 +k 2 + k 2 +p 2 +p 2^n^ k2+h 2 + k 2 Z B{k,k z ). ( B3 9) 
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Because of the scale invariance and in-plane rotational symmetry of the system, we can assume the form of z B (p 1 p z ) 
to be 

z B (p, Pz )=(^y +1 s(^), (B40) 

where A is a momentum cutoff and 5(|p z |/|p|) is an unknown function. The renormalization factor becomes Z\ oc 
A -7-1 and thus 7 + 1 is the anomalous dimension of the composite operator <fixp B . (Here 7 and 7; appearing later are 
defined so that they coincide with the definition of the scaling exponents used in Ref. [29].) Once 7 is determined, 
the scaling dimension of the renormalized composite operator 0^ B ° B is given by 

A^=°J = A + A^ B + 7 + 1 = 4 + 7. (B41) 
Substituting the form (B40) into Eq. (B39), we obtain 

| p| WNA_ A ,+i ^rn B f A dkdk z 1 Ifc| 7+1 S (M] m2) 

' Pl \\P\J m AB J (27T)3 ^(k + pY + V + k 2 + P 2 + p 2 ^ k 2 + fe 2 + fcj \\k\ ) ' { > 

In order for the integral to be infrared finite, Re (7) > —3 is necessary. Also, in order to be able to take the limit 
A — > 00, Re (7) < 1 is required. In the infinite cutoff limit A — > 00, 7 satisfies the following integral equation: 

h n + i s (M] = *™b r ^ 1 ifci 7+1 s(M\ 

' Pl \\P\J m AB J_ x (27T)3 M(fe + p )2 + fe 2 + k 2 +p 2 + p 2 Jrr^tf + + k * { |fe| J " 

Here the integral is understood to be evaluated where it is convergent —3 < Re (7) < — 1 and then analytically 
continued to an arbitrary value of 7. 

Similarly, for general orbital angular momentum I, we consider the following three-body composite operator: 

l 

0% B (x) = ZX 1 ^^ 9 * +id y ) j <f>(x)(d x +id y ) l - j i, B (x). (B44) 
3=0 

In order for 0^ BB to be a primary operator ([Ki, 0^ BB ] = 0), the coefficients Cj have to be chosen so that 

' / M Y f , m B 

y c.a I n A 

3=0 



C *\P + T7T 1 [-P+ 177 1 « P ( B45 ) 

M + mg / V M + m B 1 



being independent of the momentum q conjugate to the center-of-mass motion. In the important case of I — 1, we 
easily find c\ — — I jf-c . If we denote the anomalous dimension of such a composite operator as 7/ + 1 — I, it is 
straightforward to show that 7/ satisfies 

|„|7.+1 C (M) 4 ™ B [°° ^1 C0S (%p) Ifc| 7 ' +1 c, (\M\ fR 4 fil 

lPl l \\p\J m AB J_ x (2nf ^(k + pr + k^ + k^+p^+p^ ^k^TWTkJ l \\k\)- 1 ' 

Rescalings of the variables k — > -\f^fk and p — > \f^jfv and redefinition of the unknown function Sz(bzl/|p|) lead 
to the result shown in Ref. [29] . The scaling dimension of the renormalized composite operator 0^ BB is given by 

A% B = A* + A* B + J + (<n + 1 - J) = 4 + 71. (B47) 

By solving the integral equation (B46) numerically, we find that the anomalous dimension 71 in the p-wave channel 
I = 1 decreases with decreasing the mass ratio rriA/m B and eventually becomes complex as 71 = — 2 ± i Im (71) when 
mA/m B < 0.0351287 [29]. This implies the existence of the Efimov bound states in the p- wave ABB three-body system 
[see discussions about Eqs. (B34) and (B35)]. Furthermore, in the range of mass ratio 0.0351287 < u < 0.0660841 [29], 
the anomalous dimension is —2 < 71 < —1 and thus the scaling dimension of the three-body composite operator O^g 1 ], 
becomes 2 < A^,^ < 3. As a consequence, an additional ABB three-body resonance can be introduced and the 
resulting system provides a novel nonrelativistic defect CFT describing two-species fermions with both two-body (AB) 
and three-body (ABB) resonances in the 2D-3D mixture [see discussions about Eq. (B36)]. 

It would be difficult to determine scaling dimensions of composite operators with particles more than three. However, 
it is possible to estimate them by numerically solving the energy eigenvalue problems of H osc with the help of the 
operator-state correspondence (B16) or by using the analytic e expansions around the special dimensions d B — > 4 and 
d B -> 2 (see Sec. HIE) [52]. 
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FIG. 10: Feynman diagram evaluated in Eqs. (B49) and (B56). The solid lines are the propagators of tp A and i/jb fields and 
the shaded bulbs represent the two-particle scattering amplitude iA(pA +Pb)- 



3. Operator product expansions 

Here we consider an arbitrary effective scattering length — oo < < oo and study operator product expansions 
(OPEs) in our defect quantum field theory (1). We first work on the following OPE: 

^(*-f)^(x + |) =Y J W A<n {y)O n {x). (B48) 

n 

Here WA,n{y) arc Wilson coefficients and O n (x) are renormalized defect operators. We remind that x and y arc 
two-dimensional coordinates on the defect. We can determine the lowest three Wa.u and O n by evaluating the matrix 
elements of the both sides of Eq. (B48) between two-particle states {pa,Pb\ and \pa,Pb)- According to Ref. [57], we 
shall consider the Feynman diagram depicted in Fig. 10 that has nonanalyticity at \y\ = 0. By denoting the total 
energy and momentum as p = pa + PBi the matrix element of the left-hand side of Eq. (B48) becomes 

U\ (x - |) V* (x + f ) ) = M(p)] 2 / §^ e*» iG A{k)iG A{k)iG b{p k) 

-' 2 7 /fig: io J (2tt)3 (B4g) 



m B rriAB 



—2m B £ v 



Jl(j)") 2 e i ~M~ p " y ~\ y \v / ~ 2mAB£p 



Here A(p) is the two-particle scattering amplitude given in Eq. (7) and we introduced a shorthand notation £ p = 
Po — j^j + i0 + . If we expand the exponential in terms of \y\, the terms with odd powers of \y\ are nonanalytic at 
Ivl =0: 

/t ( y\ i ( , y\\ m B m A B ,, , 2 i^w m ABy/m B m AB ., , 2 i i , r\f 2n ^ rnN 
A\* - 2 J H X+ 2 )/ fi8 : 10 = 2^-2m B £ p A(p) & 2~, A{p) |y| + ° {y } ' (B50) 

The first term expanded in powers of y can be easily identified with the Taylor series of the left-hand side: 

-.Aipfe^-y = {^a{x)W w + J ■ (ViV^(x)) fig: i + • • • . (B51) 



m B m A B 



27ry— 2m B £ p 2 

Below we will show that the second term in Eq. (B50) can be identified with the matrix element of the defect operator; 
il) A i() B ipBipA(x) = ^ a {x)4> b {x;Q)'iI)b{x,0)'iI}a{x). 



The matrix element of ip^ipgtpstpAix) between the same two-particle states (j>a,Pb\ and \pa,Pb) is evaluated as 



/J „/,T 



{ip\ip B ip B ip A {x)) 



I - ,A( P ) l^^iG A {p-k)iG B {k) 



(B52) 



where we used Eq. (4). By comparing Eq. (B52) with the second term in Eq. (B50), we find the OPE of 
itf A (x- f) i/> A (x + Z) to be 



tf A (x - |) <Pa(x + I) = ^ A (x) + I • V^VVm(x) - mAB ^ BmAB \V\ 9o 2 ^WbMa(x) + 0(y 2 ). (B53) 

Here ffo V'^V's^ 'b^ : a( x ) is the renormalized defect operator having finite matrix elements. This result is a general- 
ization of the OPE studied in the usual 3D case in Ref. [57] to our 2D-3D mixture. In particular, the existence of 
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the nonanalytic term in \y\ implies that the two-dimensional momentum distribution of A atoms has the following 
large-momentum tail: 

PA (\k\) = f dye- ik y UUx f ) 4>Jx + f )) 

J , _ 27 V 2//any (B54) 

m A By/m B m A B , 2/ t ,t / / r \\ m \ \ 

-> ^3 (g '4>A' l pB^Bi>A{x)} a ny (\k\ -» OO). 

Here the expectation value can be taken with any state in the system, for example, at finite densities of A and B 
atoms and at finite temperature. The quantity in the right-hand side is called the contact density and given by 
(gQ r ip A tpg'ipBipA(x)} a ny — * 47T 5 ' a 2 s n A n b / ' (m ab) in the weak coupling limit a e s — ► —0. The coefficient of the large- 
momentum tail has played an important role in the exact analysis of the unitary Fermi gas in pure 3D [58] . It is an 
important future problem to investigate exact relationships in our 2D-3D mixed dimensions. 

The following OPE will be more interesting because it involves the z-direction perpendicular to the 2D defect: 

^ B (x-^,z)i) B (x + \,z)=Y,WB,n{y,z)On{x). (B55) 

n 

Wilson coefficients WB, n {y, z) and renormalized defect operators O n (x) can be determined by evaluating the matrix 
elements of the both sides of Eq. (B55) between the two-particle states (pa,Pb\ and \pa,Pb)- We shall consider the 
Feynman diagram depicted in Fig. 10 again that has nonanalyticity at \y\ = \z\ = 0. The matrix element of the 
left-hand side of Eq. (B55) becomes 



(x - |, z) ^ B (x + |, z) ) fig iq = [iA{p)] 2 J e lk y iG B (k; z)iG B (k; -z)iG A (p - k) 



m B mABA{p) e 



>- P y fJh 

J (27T 



(B56) 



) 2 k 2 - 2itiab£ p 



When \z\ ^ is fixed, the right-hand side is analytic in terms of y and therefore the OPE oiip B (x — ^, z) i\)b{x + f , z) 
is simply given by its Taylor series in powers of y. This is natural because there is no interaction with tpA(x) away 
from the 2D defect located at z = 0. 

We now set y = and study the OPE of tp B ipB{x, z) = ip B (x, z)^b{x, z) as a function of the distance from the 2D 
defect \z\ (termed defect operator product expansion). Performing the integration over k in Eq. (B56) with y = 0, 
we obtain 

(^ fl (a,z))fig:io = - mB ™ AB A(p) 2 Vi(-2\zW-2m B £ p ) 

277 ; ' (B57) 

A(p) 2 ln[2e^\z\^-2m B £ P ) +0(\z\), 



m B m A B 



2ir 

where Ei(x) = J^^e - '. We can identify the lowest order term in the right-hand side with 

_m^^^ )2 m( 2e7ENV _ 2 ^ (B58) 

where A is an arbitrary momentum scale. Therefore we find the defect OPE of ifi B tp B (x, z) to be 

4>Ub{x,z) = ^Ub(x,0)^ - H^£ln(2e^|z|A) 9 ^WbMa{x) + 0{\z\). (B59) 

Because i/j b ^b(x, z) is the density operator of B atoms, the above result suggests that the density of B atoms diverges 
logarithmically toward the 2D defect \z\ — ► 0: 

nB(\z\) = (^ b i/j b (x, z)) any -» - mB ™ AB i n \ z \ { g 2 ^ A ^ B ^ B ^ A { x )) &ny . (B60) 

The coefficient of the divergence is given by the contact density up to the mass-dependent factor. Further analysis to 
elucidate this aspect will be worthwhile. 
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4. Conclusion 

Two-species fermions in the 2D-3D mixed dimensions in the unitarity limit can be regarded as a nonrelativistic defect 
CFT. We derived the reduced Schrodinger algebra and the operator-state correspondence in general nonrelativistic 
defect CFTs. We also studied scaling dimensions of few-body composite operators and operator product expansions in 
our 2D-3D mixture. In particular, for the stability of the many-body system near the unitarity limit, we showed that 
the mass ratio has to be in the range 0.0351287 < m 2 D/"^3D < 6.35111 to avoid the Efimov effect [29]. Finally, we 
emphasize that all field-theoretical methods presented here to determine scaling dimensions and critical mass ratios 
are widely applicable to both fermionic and bosonic systems and also in the 1D-3D mixture [29] and in the usual 3D 
case [47, 52, 59]. 
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